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Abstract 



Recently, in [lj QCD sum rules for chiral partners in the open-charm meson sector have been 
presented at nonzero baryon net density or temperature up to and including mass dimension 5. 
' Referring to this, details concerning the cancellation of infrared divergences are presented and 

important technical and conceptional ingredients for an incorporation of four-quark condensates 
beyond factorization and of other mass dimension 6 condensates are collected. 

o , , 

1 Introduction 

A natural way to investigate dynamical chiral symmetry breaking is via considering the spectral differ- 
>^ ■ ence of chiral partners. In [I], chiral partner sum rules for heavy-light quark mesons in the spin and 
spin 1 channels have been presented. For heavy-light quark mesons, a crucial step in the evaluation of 
the operator product expansion is the cancellation of infrared divergences in virtue of the introduction 
of non-normal ordered condensates. Thereby, the absorption of infrared divergent Wilson coefficients 
into the condensates at non-zero densities or temperatures requires additional renormalization relations. 
In [2] , in- medium relations which are consistent to the vacuum case have been given. However, only the 
pseudoscalar case has been investigated there and, moreover, a general proof to all orders in the mass 
dimension is still missing. Chiral partner sum rules may, therefore, serve as a consistency check of the 
relations derived in [2]. 

In view of the dynamical chiral symmetry breaking of quantum chromodynamics and order parame- 
ters of the symmetry restoration (phase) transition, QCD sum rules for heavy-light quark pseudoscalar 
mesons provide the amplification of the numerical impact of the chiral condensate in virtue of the heavy 
quark mass. Hence, four-quark condensates, which are of mass dimension 6, are numerically not as im- 
portant as in the case of mesons consisting of light quarks only. However, in order to determine their 
numerical impact one has to take them into account. The evaluation of Wilson coefficients for four-quark 
condensates in case of heavy-light quark sum rules is an elaborate task. A common approximation is the 
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vacuum saturation hypothesis, which leads to a factorization of four-quark condensates into squares of 
the chiral condensate. Once applied, the factorization significantly simplifies the calculation. However, 
as only certain combinations of four-quark condensates may serve as order parameters of chiral symme- 
try the factorization of all four-quark condensates in terms of the chiral condensate is not necessarily 
consistently Therefore, the extension of mass dimension 5 in-medium heavy-light quark QCD sum rules 
to mass dimension 6 requires the consistent evaluation of four-quark condensate contributions. With 
respect of this task, chiral partner sum rules provide assistance as only chirally odd condensates do not 
cancel out and the number of terms which have to be calculated is significantly reduced. 

2 Chiral partner heavy- light quark meson correlators 

We consider the currents 

j S 0) := ?i(ar) <h(x) , f(x) := q^x) 175 q 2 (x) , (la) 

jj( x ) Qi( x ) 7/i Q2(x) , j A (x) := qi{x)^ fl q 2 (x) (lb) 

and the corresponding causal correlators 

n( s < p )(g) = i J d 4 a;e^(T [j (s ' P) (^)j (S ' P)t (0)] ) , (2a) 

U%> A \q) = 1 J d 4 xe^(T [j^ A \x)j^ (0)} ) , (2b) 

where T[. . .] denotes time-ordering and (. . .) means Gibbs averaging [3J. In the rest frame of the 
nuclear medium, i.e. n = (1,0), and for mesons at rest, i.e. q = (1,0), the (axial-) vector correlator can 
be decomposed as 

= (^f - 9 S) 4 V ' A) (g) + q ~f^ A \q) (3) 

with U.^' A \q) = 1 g^^j H$,' A \q) and U^' A \q) = ^q^q v ^-% ,K \q) written manifestly covariant. 

The transversal projection for a vanishing light quark mass is given by 



4" A (g) = -0n p - s (g) - V- A ( g ) - f^iMi) , (4) 

with n( v > A )(g) = g^U^' A \q). 



3 Cancellation of IR divergences 

In lowest-order perturbation theory the operator product expansion within the background field method 
[H [5] for the difference of chiral partner correlators reads 



n p - s (°)(g) = -i j ^{Tr^ip + q^Tr^p)] 



+ ^Ti D [S 2 (p + q^TrviSiip)^] +Tr D [S 2 (p + ^Tr^phs]} :> , (5a) 



1 Note, that the factorization of four-quark condensates into squares of the chiral condensate is exact in the limit of 
infinitely many colors. 



2 



n ' r 

+ g 2 ra 5 « (Trotr^^-g)]) ^ n g 2 :) (5b) 

n V - A(0) (<7)=* | 2Tr c {Tr D [5 2 (p + g)]Tr D [5 1 (p)] 

-Tr D [5 2 (p + g)75]Tr D [5 1 (p) 7 5]} :> , (5c) 

nV_A(2) (?) = -E^f 2 E (:?iKr^(Tr D [r5 2 (g)]) gi 
n ' r 

+ g 2 r^" (Tr D [r5!(-g)]) D 5n g 2 :) , (5d) 

where Ttc,d denotes trace w.r.t. color or Dirac indices. D^(x) = d^ — igg/^x) is the covariant derivative 
and an arrow indicates whether it acts to the left or right. For sake of a concise notation we have 
defined D$ n = D ai . . . D an (with an analog notation for the partial derivative). The quark propagator 
is denoted by iSij{p) = J d A xe ipx (fi|T [qi(x)q~j(0)] \ fl), \fl) is the ground state of the strong interaction, 
0>k = ^[7^,7^/2 and : ... : means normal ordering w.r.t. the perturbative vacuum |0). Eq. (J5]) 
stems from the application of Wicks theorem to ([2]) in lowest order of the perturbative expansion, 
n = n < -°' ) + n^, where corresponds to the term where the quark fields of the currents ([1]) are 
contracted to propagators and to the term with 2 not contracted quark fields. There is no flow of 
hard momentum in the term with four not contracted quark fields. Hence, it does not contribute. 

In P, it has been proven that TrD[«Si(p)r] = for mi = 0, p 2 7^ and T G {1, c^, 75}. Hence, a 
superficial view on ([5]) may tempt to the conclusion that II is zero for chiral partner operator product 
expansions of heavy-light quark meson currents. On the other hand it is clear that only matrix elements 
of chirally odd operators may enter the OPE, whereas gluon condensates are chirally even. In this sense, 
the cancellation of n^ ) is in line with naive expectations. But the introduction of non-normal ordered 
condensates via n^ 2 ^ would introduce gluon condensates together with infrared divergences. Taking this 
as a heuristic argument for non-zero n(°) raises the question of the precise cancellation of these terms 
and terms added by introducing non-normal ordered condensates. Indeed, from the in-medium OPE of 
D mesons up to mass dimension 5 [21 E] it is known that the medium specific divergences are canceled 
because of the renormalization of (: d^^Dyd :). But this term is chirally even (formally, it is the matrix 
element of a vector current) and does not enter §5§. Moreover, it is known that introducing non- normal 
ordered condensates, in order to cancel infrared divergent Wilson coefficients of gluon condensates, 
leads to additional finite gluon contributions. Clearly, these have to cancel out in case of chiral partner 
sum rules. Again, this can be taken as a heuristic argument that only those mass divergences can 
remain m n(°) which are cancelled by chirally odd condensates. Two questions have to be answered. 
Do all infrared divergences cancel out? And does the renormalization procedure introduce chirally even 
condensates? Thus, a careful analysis is mandatory to prove that the obtained results are infrared 
stable and that the renormalization procedure is consistent. 

The quark propagator in fixed-point gauge (x — xo) M ^ M (a;) = [5] can be written as S(p) = 

J2n=o Sin) (p) with S^{p) = (-l)S^(p) (7 A) S^(p) = (-l)S®(p) (7/) S^(p) . A denotes a 
derivative operator which arises due to the Fourier transform of the perturbation series for the quark 
propagator in coordinate space from the gluonic background field ^/^ A^ = Xl^Lo ^jf 1 with Ay = 
~ { -^W (jDai ' ' ' dU ° ai ■ ■ ■ 9an » where % v = i[D^ D v ) /g = G%t A is the gluon field strength 

tensor and g = \/4na s is the coupling. t A are the generators of the color group and A = 1, . . . , N 2 — 1. 
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For the contributions to fl5a|) up to mass dimension 5, which restricts the quark propagator to next-to- 
next-to leading order and the gluon field to lowest order, the following traces are evaluated 

Tr D [S^(p)] = -^- 2 , Tr D [S^(p) l5 a, u ] = , Tr D [S^(p) l5 ] = , (6a) 

Tr D [S^(p)j 5 a^] = -ig *™ e, UKX ^ x , Tr D [S^(p)] = , Tr D [S (1 \p) l5 ] = , (6b) 

\P Tfl ) 

Tr D [^)( P )]=8/^^^n. (6c) 

It is not necessary to consider traces of the second order quark propagator with 75 and a^ u . Up to 
mass dimension 5 they can only be multiplied with their lowest order counterparts fl6aj) . which are 
zero. A combination of second-order and first-order propagator leads to mass dimension 6 terms. Using 
TYd [Si 7507^] TrD [S^sO"^] = Tro [-Sicr^] TrD [SW^] the following contributions have to be considered 

n p - s (°)( g ) 

^jj<: ^Trc[Tr D [S(°)(p + g)]Tr D [Sf (p)] + Tr D [Sf (p + 9 )]Tr D [sf (p)] 

+Tr D [Sf (p + g)]Tr D [Sf (p)] + \ttb[S^\ P + q)a, u ]Tr D [S d 1 \p)a^ :) . (7) 

Due to Tro [5^75] = Tro^ 1 -^] = for arbitrary quark masses, there is no 75 contribution up to 
this mass dimension. Likewise Tro [S^^cr^] = Ttd^W] = 0. Lorentz invariance requires that terms 
with only one gluon field are zero. Therefore, a first-order quark propagator must be combined with a 
propagator of at least the same order. Keeping both quark-masses m rf c finite these four terms give rise 
to the following integrals 



n p - s (°)( g ) 

_. r d 4 p 1 4m c Am d 
% J (2tt)^' 2{p + q) 2 -mlp 2 -m 2 d T ° [ Cj 

+ ~Q~ T7 Z ^2 2U~2 2P P G ^ G a Tr C [t t 

2 [{p + q) A — 171^1 p — m d p L J 

, 9 2 4m c 8m d ^, ar ,A r B !/ T \+a+b^ 

+ -2{p + q y-ml[p>-mlf P ° Tr ° V 1 J 



+ J [{p + q) 2„ m 2 c] 2 [p 2_ m 2 ? G ^ G T ^C [t t \:}. (8) 



Analyzing the integrals in euclidean space in terms of the integral 

o 2 2\ f 1 — a V 

hjk(q ,m d ,m c ) = da— TvZ^T ( 9 ) 



reveals the following results. The first term is cx mrfm c I 0i o,o(9 2 ) w^, m^). The second term generates two 
terms which are cx m<im c Io,3,2((7 2 , wi^, raj!) an d oc mdm c l2,3 t 3(q 2 ,m d ,ml), respectively. The last term is 
oc mdm c li t i^{q 2 ,m d ,ml). These terms are all zero in the limit m d — > 0. On the other hand, the third 
term does not vanish for m d — » 0. It causes a term oc mdm c l 3 fi : 2(q 2 , m d , ml) , which diverges with m^ 1 . 
Projection of the Lorentz indices leads to 

<: :) = G 2 :) - I (g m - 4^) (: - :) , (10) 



where i> M is the medium four velocity, (vG) 2 = v^v v g a ^G A a G A p and G 2 = G^ v G Av ^ . Note that the 
diagonal elements of (fTUj) are vacuum specific, whereas the medium specific contribution is traceless. 
The integral of the third term of (j7|) can be evaluated in euclidean space: 

d^p m d m c (p - q)^(p ~ q) a 



(27r) 4 [(p — q) 2 — m 2 ,] 4 p 2 — m 2 

^ \ (^ l3 > - 2 ^ 2 ' m ^ + 2 ^^ I 3,2,3(g 2 , ml, m*)) , (11) 

where I 5A3 - 2I 4)0) 3 + 13,0,3 = 13,2,3 has been used. By virtue of l ijk (q 2 ,m 2 d ,m 2 c ) -> l jik (q 2 , m% m 2 ) the 
corresponding integral of the second term in ([8]) can be derived. The limit m d — > for both terms in 
(HU) is 

lim m d l 3A2 (q 2 ,m 2 d ,m 2 c ) = — — — ^ — - , (12a) 
md^o m d q 2 + m~ 

lim m d I 3i 2, 3 (g 2 , rrid, m 2 ) = , (12b) 

m d ->0 

where (I12al) reveals the famous infrared singularity from vacuum D meson sum rules. As the medium 
specific contribution to ( ITU]) is traceless with respect to Lorentz indices and due to the vanishing of the 
second term in (II ip there is no medium specific infrared divergent term. Hence, the only terms that have 
to be absorbed into the condensates by virtue of the introduction of non-normal ordered condensates 
are vacuum specific. This is in line with the cancellation of (: dj^D^d :) in (jSJ) which would have to 
absorb the medium specific divergences. Moreover, infrared divergent terms which enter through the 
medium specific gluon condensate, therefore, must enter the 7 M or 757^ parts of the quark propagator. 
Owing to (ITDjl . (ITT]) and (|T2|) the limit m d — > of ([7]) in euclidean space is 

n p - s (°)( g ) 

= ~/ (0I<: Tr D [^)(p + g)]Tr c , D [5f(p)] :) ^ -^(:^G 2 :)I 3 , ,2(g 2 , m% m 2 c )m c m d 

= -^<:-G 2 :)^^, (13) 
o 7r m d q 1 + m~ 

where it has been used that the lowest order quark propagator is a unit in color space and the additional 
imaginary unit stems from the Wick rotation (W.R.). 

The normal ordered chiral condensate in lowest order of the light quark mass enters (I5bl) via the 
following expression [1] 

^ (2 \q) = l(-.dd:)TT D [S c (q)} . (14) 

Here, Tr^ [S^oO] is the Wilson coefficient of the chiral condensate. Expressing the normal ordered 
condensate by the non-normal ordered condensate 

(: dd :) = (dd) Tr c , D [S d (p)} :) (15) 



leads to 



n S) (2) (<?) = \(dd)Tr D [S c (q)\ + ~Tr D [S c (q)\ * J ^0 Tr c , D [S d (p)] :) . (16) 
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Figure 1: Soft gluon momenta four quark condensate contributions from higher order quark propagators 
and quark field expansion. 

Despite the striking similarity, revealed by the projection onto elements of the Clifford basis in (jSJ), of 
the third term in ([7]) which has to be canceled and the additional term in (fl6|) . a general proof cannot 
be given and the necessity for a precise evaluation is obvious. Evaluation of (fl5|) up to order a\ in 
MS-scheme gives 



with the renormalization scale /i. Insertion into ( fT6l) and in virtue of the limit rrid — > the anticipated 
cancellation of infrared divergences in terms of light quark masses is revealed 



Adding ( Tl2l) and (fl8l) the infrared divergent term cancels out. Furthermore, Eq. ([5]) and the explicit 
evaluation of the renormalization of normal ordered condensates [2] shows, that in chiral partner sum 
rules up to and including mass dimension 5 in the limit — > only the chiral condensate mixes with 
other condensates by virtue of introducing non-normal ordered condensates. Thus, apart from the term 
which cancels the infrared divergence no additional chirally even terms enter and the final operator 
product expansion is chirally odd. 

So far, the investigation was carried out for the spin case. Fortunately, as the only formal difference 
between the spin and spin 1 case is the cancellation of the projections of quark propagators onto the 
tensor a^ u in the spin 1 case and a global factor —4, see ([5]), from the previous evaluation it is clear 
that the terms of interest are the same in both cases up to mass dimension 5. 

4 Comments about four-quark condensates 

We now turn our attention to four-quark contributions in fl2]). The typical four-quark condensates for 
the a, corrections are of the form 



where Y e {1, 7 M , cr^, 757^, 75} is a Dirac matrix and r G U(iVf) ® U(iV c ) is a matrix in flavor and 
color space (respectively, N c $ is the number of flavors and colors). L denotes an n-fold tensor in 
Minkowski space constructed by suitable combinations or contractions of {v^, g^ u , e^ap} such that ( Fl9l) is 
a Lorentz scalar, in particular even w.r.t. parity. t^ ua p is the Levi-Civita symbol, a pseudotensor. Their 
diagrammatical representation can be found in Fig. [1] for the contributions with soft gluon momenta 
and in Fig. [2] for contributions of hard gluon momenta. 

These can be grouped into such ones which are already present in the vacuum sum rule and others 
which are specific for the medium. For the medium specific condensates the medium velocity v enters L. 




(17) 




(18) 



(^rr^rVV)L , 



(19) 
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(a) 



(b) 



Figure 2: Hard gluon momenta four quark condensate contributions from higher orders of the pertur- 
bative expansion. 



It should be emphasized that four-quark condensates (cf. [3, El El HOI EH] for a complete classification) 
enter in different combinations the light vector meson (cf. PHHE]) or the nucleon sum rule (cf. [TUf [7]). 

Soft gluon momentum contributions arise from non-local quark fields or higher order quark prop- 
agators. Expanding non-local quark fields up to third order in the covariant derivative and using the 
equations of motion for the gluon field strength tensor gives rise to the diagrams depicted in Figs. 



1(a) and 1(b) Inserting quark propagators with a next-to- leading order gluon line attached and using 



the equations of motion for the gluon field strength tensor leads to the diagram shown in Fig. 1(c) 
Soft gluon contributions in chiral partner sum rules are proportional to the medium velocity v and are 
therefore absent in vacuum. 

Contributions with hard gluons arise from the 0(a s ) perturbative correction. There are vacuum as 
well as medium condensates. 



5 Summary 

In summary, we try to extend the chiral partner QCD sum rules derived in |Tj in zero order of a s . Higher 
order of ot s are related to higher-dimensional condensates. Most notable the four-quark condensates, 
which are of mass dimension 6, enter such difference sum rules. Abandoning the vacuum saturation 
hypothesis introduces a large number of different vacuum and medium four-quark condensates, but it 
is necessary to keep a clean separation of chirally odd and even terms. 

Furthermore, gluon condensates of the form (G 3 ) are of mass dimension 6 as well. Possible in- 
frared divergences linked to them and their cancellation for finite temperatures or densities have to be 
investigated carefully. 
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